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1. Introduction 

A variety of well established methods can be applied to investigate diffusion in solid 
state systems [1]. However, studies of diffusion mechanisms on the atomic scale are still 
challenging. Atomic scale X-ray Photon Correlation Spectroscopy (aXPCS) can provide 
new insight into this field. 

XPCS is the X-ray analogue of the very successful technique of dynamic light 
scattering or Photon Correlation Spectroscopy [2]. PCS is well established to 
characterize the equilibrium dynamics of soft condensed matter by determining the 
intensity autocorrelation function of the scattered laser light as a function of delay time 
and wave- vector transfer. Similarly, XPCS |3J, using coherent synchrotron radiation 
instead of visible laser light, is being widely used for studies of soft matter dynamics 
where objects are in the nanometre range (see e.g. [U El EJ [7] ) . Attention has also been 
directed towards studies of metal systems as shown by the feasibility test on the (100) 
Bragg peak of CU3AU [8] or of critical dynamics in FesAl [9]. Because of the short 
X-ray wavelength, XPCS seems to be predestined for studies at atomic length scales. 
These types of experiments have, however, never been reported in the literature until a 
successful demonstration of the capability of XPCS to resolve single atomic motion in 
condensed matter [10] . Particularly, we have clarified the mechanism of atomic dynamics 
in CugoAuio [El E] and studied dynamics in a metallic glass [12] . 

The goal of this paper is to describe in detail the atomic motion of dilute platinum 
atoms in a nickel matrix at relatively low temperatures. The Nig7Pt3 intermetallic alloy 
is a perfect candidate for studying diffusion of small quantities of dilute atoms due to 
its outstanding intensity of diffuse scattering even for low solute concentration levels. 
Moreover, some profound questions can be addressed regarding interactions between 
solute atoms and interactions of solute atoms and vacancies. It is well known that 
vacancies are vehicles of diffusion in alloys [1] , but the interactions between vacancy and 
solute atom can radically influence the motion of the latter. Our objective is to analyse 
these interactions with the help of aXPCS. These questions were addressed in the past by 
a number of methods like Mossbauer spectroscopy, nuclear resonant scattering or quasi- 
elastic neutron scattering [13]. The elementary diffusion steps were investigated for self 
diffusion [14] . for impurity diffusion with consideration to the interaction with vacancies 
[15] , for hydrogen in intermetallic hydrides [16] and in intermetallic B2 [T7 | IT8 | IT9] and 
DO3 [20J phases. All methods mentioned above have in common that they are limited to 
only few selected isotopes and have a poor energy resolution which forces measurements 
at very high temperatures (apart from hydrogen). The time-domain technique of aXPCS 
is, however, free from these restrictions. 

The important aspect of this paper is studying diffusion in single crystal and 
polycrystalline samples of the same composition. For methods studying spatially 
resolved atomic motion like aXPCS, single crystal samples provide the possibility to 
investigate different crystallographic directions and therefore more information. It is, 
however, sometimes difficult or even impossible to obtain single crystals of satisfactory 
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quality. Showing that studies on polycrystals can, in whole or in part, provide 
information on a comparable level when using a coherent method, a new field of materials 
could be investigated with aXPCS. Ni 97 Pt3 seems to be a predestined object for such a 
study. 

2. Theory 

Atomic structures can be investigated at their elemental scale. Still, monitoring a single 
atom in a bulk sample on time scales of atomic diffusion it is not possible in real space. 
It is however possible to measure dynamics in reciprocal space which corresponds to 
these dynamics in real space. 

Usually, reconstruction of real space from a scattering image is very challenging 
due to the phase problem. Fortunately, one can ignore the details of the exact structure 
factor, but use the time fluctuations of the intensity pattern to study dynamics of the 
material with XPCS. This is possible because the intensity autocorrelation function 
g^ 2 \q,t) and the autocorrelation function of the electric field g^\q,t) are connected by 
the Siegert relation for Gaussian scattering processes [21]. Therefore it is sufficient to 
measure the intensity autocorrelation function to gain information about the dynamics 
in real space 



where I(q, t) is the observed intensity at wave-vector q and time t and the brackets (...) 
denote an average over the ensemble. In XPCS and in many other types of experiments, 
the time average of the investigated quantity is calculated from the measured data. 

XPCS is the time Fourier-transformed counterpart of coherent quasi-elastic neutron 
scattering (QENS). It can be shown that the correlation time r(q) corresponds to the 
energy line broadening [22] T(q) via r(q) oc T(q)^ 1 . Therefore concepts developed for 
QENS can be adopted to XPCS. 

For a coherent method like aXPCS two atoms of the same species exchanging their 
positions do not affect in any way the scattering pattern, as the configuration of the 
coherently illuminated volume before and after the jump is exactly the same. Compared 
to non-coherent methods it is, however, necessary to account for the short-range order, 
i.e. for the influence of the surroundings on the movement of each distinct atom. This 
phenomenon has first been described qualitatively for coherent neutron scattering in 
liquids by de Gennes [23]. The physical principle for the so-called de Gennes narrowing 
is based on the fact that peaks in the structure factor of liquids occur at wave-vectors 
corresponding to the most probable interatomic separation, i.e. they are responsible for 
the highly correlated and long living atomic arrangements. The same argument holds 
for the short-range order intensity (Jsro) in alloys where the higher diffuse intensity 
between Bragg peaks corresponds to a longer lifetime of certain configurations of solute 
atoms. 
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The simplest case of Lorentzian line-shapes with energy-domain methods 
corresponds to an exponential decay of the intensity autocorrelation function in time 

g {2) (q,At) = l + /3exp(-2At/<r(g)). (2) 

The intensity autocorrelation function for perfectly coherent radiation would decline 
from a value of 2 to 1. However, because of an only partially coherent X-ray beam 
in the experimental setup the amplitude /3 is diminished. The correlation time r(q) is 
given by 

En P n Y,Aa nl i 1 ~ eX P ( J Q ' Aa nj)) 

This follows from the theory describing the short-range order and coherent dynamics in 
a linear approximation as derived by Sinha and Ross |25j . It is discussed in more detail 
in Leitner and Vogl [11]. tq is the average time between jumps of an atom and therefore 
corresponds to the inverse of the average jump frequency v^ x . The probability for an 
atom to jump to one of the lattice sites of a certain neighbouring shell is given by P n , 
with n being the number of the appropriate coordinate shell (1 for nearest-neighbour, 2 
for next-nearest-neighbour and so on). The vector a n j points from the origin to the jth 
site in the nth shell. It can be easily seen that the inverse correlation time is the sum 
over the inverse correlation times for each shell corresponding to the particular jump 
mechanism. Therefore it can also be written as r{q)~ l = En r «(^) _1 - 

In a polycrystalline sample one has to average over all orientations of q with respect 
to the crystal lattice 

W(q, At) = 1 + A / dq • exp (-2At/r(q)) . (4) 

It is essential for the method at hand to be able to fit the data set with a function 
depending on At. As there is no analytical function which satisfies equation (j3J), 
an approximation has to be made. A similar approach was introduced by Chudley 
and Elliott (26] for liquids. Exchanging integration and exponentiation, an average 
correlation time r(q) can be defined as r^) -1 := (47r) _1 f\* =q dq ■ r(g) _1 . For this 
expression an analytical solution can be found (Isro(q) was a l so integrated over \q\ = q): 

-I \ lsRo{q) fK \ 

T{q) = T 7 r-. (5) 

p 11 — sm (g- A "rcj) 
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Using equation (jSJ) the intensity autocorrelation function can be written as: 

g {2 \q,At) = l + /3exp(-2At/r(g)). (6) 

Numerical calculations show that equation (jBJ) is a very good approximation for 
equation (0J away from Bragg peaks. 
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Figure 1. Autocorrelation function g^ 2 ' (q, At) in Ni-Pt single crystal orientated with 
(110) direction parallel to the incoming X-ray beam with azimuthal angle <f> = 33.2° 
at 855 K (upper plot) and at 815 K (lower plot) fitted with equation ([2]). 

3. Experimental 

Nickel and platinum, both of 99.98% purity, were melted by induction heating in a 
"cold boat" in argon atmosphere. For the preparation of the polycrystalline sample the 
resulting ingot was cut and rolled between stainless steel plates down to a thickness of 
17 |J.m. It was then pre-annealed at a temperature of 773 K for 22 h in a vacuum of 
about 10~ 7 mbar. In order to obtain sufficiently large single crystalline domains the 
rest of the ingot was heat treated at 1573 K for 140 h resulting in crystallites of about 
about 5 mm diameter. One domain was oriented, cut with a wire saw and mechanically 
polished to the final thickness of about 18 [xm. The composition of the sample was 
analysed by energy- dispersive X-ray spectroscopy (EDX) showing 3.0(5) at.% platinum 
content. 

The experiment was performed at beamline ID10A of the ESRF using a standard 
coherent setup with 8 keV photons at a monochromaticity of AE/E = 10 -4 . Rollerblade 
slits were used to select an 8 x 8 \im 2 beam section and speckle patterns were recorded 
with a direct-illumination CCD camera with 1024 x 1024 pixels at a size of 13 x 13 [im 2 . 
Prior experience has shown that the CCD camera placed close to the sample (about 
0.6 m) is advantageous to maximize the measured signal. It was mounted in a way that 
allowed it to be moved along a scattering angle 20 and an azimuthal angle 0. Note 
that in a scattering experiment the length of the scattering vector is connected to the 
scattering angle 2B via q = 2\k\ sin(2B/2), with k being the wave vector of the incoming 
beam. As one can move freely on the Ewald sphere in reciprocal space by changing the 
detector position in real space, the intensity autocorrelation function was measured for 
several g(2©,0). In the polycrystalline sample only the scattering angle 20 (i.e. \q\) 
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can influence the correlation time. Therefore measurements at arbitrary azimuthal 
angles are equivalent. For the single crystal sample an additional second scan for fixed 
\q\ and different azimuthal angles <fi was performed. Measurements were performed 
in transmission geometry in a resistively heated vacuum furnace (p < mbar). 
The temperature was stabilized by a PID temperature controller within a range of 
0.1 K. The measurement time at a fixed scattering vector q was between 20 and 40 
minutes at a frame rate of one frame every 10.9 seconds (exposure plus readout time). 
Due to the limited intensity of the synchrotron source and low scattering intensity at 
large scattering angles in the diffuse regime it is beneficial to make the data collected 
by the CCD camera subject to the so-called droplet algorithm [27] to detect single 
photon events. Our implementation of the droplet algorithm is described in [24] . These 
corrected images were evaluated according to equation (CQ), finding the correlations in 
the fluctuating speckle intensities (see figure [T|) where the autocorrelation function was 
averaged over all pixels to obtain sufficient statistical accuracy. 

All measurements were carried out in the diffuse regime, away from Bragg peaks. 
Any kind of strain broadening / asymmetry resulting from dislocations, planar faults 
etc. are constrained to q- regions very close to the Bragg reflections [28]. Therefore 
such effects do not influence the intensities obtained in an aXPCS experiment. The 
number of vacancies relative to the number of atoms in Ni at 830 K is in the order 
of 10 -8 [29j. Therefore the total scattering power of the vacancies and their induced 
displacements is negligible. An estimation using the Debye model at this temperature 
gives a contribution in the order of 10 -3 originating from thermal diffuse scattering. 
This contribution can therefore also be neglected. As a result the diffuse intensity is 
only due to the stochastic occupations of lattice sites by Ni or Pt, and the measured 
correlation decay therefore corresponds to exchanges of atoms. 

Small absolute values of reciprocal vectors correspond to very long correlation times. 
Instabilities of the experimental setup result in additional intensity variations. This can 
lead to fault when overlapping with the slowly decaying autocorrelation function caused 
by the sample itself but is not important if the time-scale of both processes is clearly 
different. Correlation times measured at very small 20 angles were therefore not used 
for the evaluation of diffusion models. 

It is vitally important for aXPCS measurements to be performed in thermal 
equilibrium. To ensure that order relaxation is negligible, multiple measurements were 
made repeatedly at different detector positions for each sample. Usually the sample 
relaxes faster when starting from a higher temperature and shows relaxation times of 
about one hour. After that time the single crystal stayed stable during the whole 
experiment. The polycrystalline sample, however seemingly relaxed again after about 
4 h at 830 K, leading to still longer correlation times. 

Note that the correlation time as measured in an aXPCS experiment quantifies 
the timescale over which the atomic arrangement changes by a sizeable amount. Even 
though effects such as enhanced grain boundary diffusion can affect the mass transport 
over macroscopic distances significantly, the contribution of the grain boundaries to the 
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scattered signal and therefore to the correlation decay is negligible. As a consequence, 
aXPCS is inherently sensitive only to bulk diffusion. 

To investigate whether there were crystallite growth effects in the polycrystalline 
sample during the experiment, the sample was studied with electron microscopy. 
Unfortunately the results were inconclusive. Therefore we assumed that such effects 
can be neglected as long as no relaxation in the correlation time is detectable. 

4. Monte Carlo simulations 

To simulate a simple diffusion model, Monte Carlo simulation of diffusion is a valuable 
approach. Two kinds of atoms were arranged in a proportion of Npt/N^ = 3/97 over 
an fcc-type lattice with periodic boundary conditions. One vacancy was introduced by 
randomly emptying a lattice site. The simulation consists in iteratively choosing one of 
the 12 lattice sites of the vacancy's first coordination shell randomly, then, according to 
the Metropolis algorithm, calculating the probability Pi^f = exp (— AE/ksT), where 
AE = Ef — Ei denote the system energies calculated with pair potentials for nearest- 
neighbours before (i) and after (/) the jump, respectively, where ks is the Boltzmann 
constant and T the absolute temperature. A decision to perform a jump or not was 
taken by comparing ft-*./ with a random number. A unit time (MC step) was defined as 
a number of trials equal to the number of lattice sites. After equilibrating the system, 
the actual diffusion simulation was started. After every half Monte Carlo step a discrete 
Fourier transform of the lattice was carried out and the amplitudes were used to calculate 
the scattering intensity and the amplitude autocorrelation function g^(q, t) for certain 
q values. Finally, using Siegerts relation, the intensity autocorrelation function was 
calculated and compared with the experimental correlation time. 

The same method was used to calculate short-range order intensities from a 
reciprocal space map of 64 3 voxels with an atom configuration of 4x32 3 atoms averaging 
the intensities over 10000 Fourier transforms. 

5. Results 

In order to distinguish between different diffusion models, measurements of intensity 
autocorrelation times in polycrystalline and/or in single crystal samples alone are 
insufficient. As shown in section [2] short-range order intensity plays a crucial role in 
the interpretation of the measured correlation times. 

5.1. Short-range order 

The particular importance of local atomic arrangements in alloys for the interpretation 
of aXPCS measurements follows from equation ([3]). Detailed data about SRO like in the 
case of CugoAuio [30] is rarely available. SRO must be actually measured for exactly the 
same alloy composition and temperature as the sample under study. It is very difficult to 
deduce short-range order data from experimentally measured intensity data. A detailed 
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description of the diffuse scattering theory can be found in several articles and books (see 
e.g. pH [32] and literature cited therein). The factor which actually makes the analysis 
of diffuse scattering so difficult and time consuming is the separation of SRO scattering 
(which is of interest to us) and atomic displacement scattering. Simulations showed 
that in a system with a very high scattering contrast and low solute atom concentration 
like Nig 7 Pt3, the magnitude of those two terms can be of the same order. Also terms 
of higher order in the expansion of diffuse intensity, e.g. caused by static displacements 
(Huang scattering) or thermal diffuse scattering, can play a role. 

Phase diagrams of the Ni-Pt system [33] show three ordered equilibrium phases 
similar to Cu-Au: NiPt (Llo structure), NisPt, and NiPt3 (both LI2 structures). There 
was a weak tendency found for a C^Au-type ordering for Ni^^-Pt^,. with x = 0.2 
composition [34] which is also indicated by the diffuse intensity for x = 0.03 as shown 
in the inset in figure El There is, however, no data for short-range order configuration 
on a system with x = 0.03 available in the literature. Note that our main goal remains 
an atomistic diffusion mechanism and the details of the diffuse scattering intensity are 
not vitally important. Therefore we choose two simple models: (ml) a model with no 
interaction between the atoms and (m2) a model with strong Pt-Pt repulsing force. 
The latter was chosen because a Ll 2 ordering type tendency resulting from strong Pt- 
Pt nearest-neighbour repulsion was assumed to be more likely than from Pt-Pt next- 
nearest-neighbour attraction or further ranging interactions 0. The two model systems 
ml and m2 were simulated using a Metropolis Monte Carlo algorithm (see section HJ). 
From the atomic configuration we calculated the Warren-Cowley parameters |35j. In 
case of 3% solute atoms with strong repulsive force the calculated parameters were 
cti = —0.0309, CK2 = 0.0042 and 0:3 = 0.0018. These parameters can be used to calculate 
the short-range order intensity in Laue units according to 



where Aa nj again is the family of vectors for a certain shell as discussed in section [2j 
As shown in figure [2] the so calculated short-range order intensity is in good agreement 
with the intensities calculated directly from the Fourier transform of the real-space 
configuration. 

5.2. Diffusion mechanism 

The main objective of this paper is to discuss the mechanism which governs diffusion in 
a Ni-Pt intermetallic alloy. An essential advantage of a microscopic (atomistic) method 
in diffusion investigation is the possibility to identify such a diffusion mechanism. In 
order to verify which diffusion mechanism operates in a specific intermetallic phase, it is 
necessary to develop a mathematical model for the diffusion process. This mathematical 
idea is generated from the knowledge of the system or it is simply a right guess. The 

X A more thorough investigation of the interaction parameters in the system Nig7Pt3 based on ab initio 
calculation will be subject to a following paper. 
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Figure 2. For a single crystal in (110) orientation with azimuthal angle <fi — 33.2° the 
short-range order intensity was calculated without interaction between the atoms (red 
triangles) and for strong Pt-Pt repulsion (blue disks). For the latter the short-range 
order was also calculated from Warren-Cowley parameters gained from the simulation 
(grey line). Inset shows measured diffuse intensity as a function of the wave- vector 
transfer. 

mathematical model can be (at least in certain cases) solved analytically and yields the 
wave- vector dependent correlation time in a single crystal sample (see equation (E])), or 
the correlation time dependent on the magnitude of the wave- vector in a polycrystalline 
sample (see equation ([5])). A model is determined by the jump probability parameters 
P n . The only free parameter fitted to the experimental data is the inverse average 
residence time Tq 1 which corresponds to the jump frequency. 

Needless to say it is not particularly difficult to predict the atomistic mechanism of 
diffusion in Ni-Pt. It is simply an exchange of a nickel or platinum atom with a vacancy. 
As the scattering length of platinum considerably exceeds that of nickel and vacancies 
are very rare in the alloy at the temperatures of measurement, the measurable speckle 
fluctuations result predominantly from the exchanges of platinum atoms with nickel 
atoms. This very simple picture enables studies of further subtleties of the diffusion 
mechanism, such as the influence of interaction between solute atoms or between solute 
atoms and vacancies. 

A simple and useful tool when discussing a diffusion model involving a vacancy 
is the so-called encounter approximation [36J. It has been successfully applied in the 
interpretation of diffusion in nuclear magnetic resonance and in Mossbauer spectroscopy 
studies [TTt IT8] . An encounter is here defined as the sum of all exchanges of one and 
the same atom with one and the same vacancy. As the number of vacancies is much 
smaller than the number of atoms, the time between one and the same vacancy jumping 
twice is much smaller than the time between consecutive jumps of one and the same 
atom. Details during an encounter are thus "invisible" for aXPCS and only a net result 
of one or more exchanges with a vacancy, which can carry an atom to nearest or to 
further neighbour shells, is accounted. One can assume that encounters are mutually 
independent and calculate e.g. diffusivity as a simple random walk of the vacancy. 
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In an fee lattice the first neighbour shell is at a distance of a/y/2, the second of a, 
the third of a^3/2 and fourth of ay/2, where a is the lattice constant. The nearest- 
neighbour shell can be reached by only one exchange between a solute atom and a 
vacancy. The second, third and fourth shells can, however all be reached after two 
exchanges with a vacancy. Still farther shells (n > 5) need three or more exchanges 
and can therefore be neglected. Also encounters that do not lead to an effective jump 
are invisible for an aXPCS experiment and are therefore neglected from here on (P 
= 0). We used this to re- normalize the jump probabilities for self diffusion given in 
[37] and gained jump probabilities of: Pi = 92.6 %, P 2 = 2.4 %, P 3 = 3.9 % and P 4 
= 0.8 % (the neglected probability for further jumps sums up to ^2 n P n >5 = 0.3 %). 
This encounter model (ml), however, describes only systems with negligible interactions 
between matrix atoms, solvent atoms and vacancies. Therefore we extend the encounter 
model according to equation ([3]), using the SRO intensity for strong Pt-Pt repulsion as 
discussed in section I5TT1 We use this as a second model for the description of Pt diffusion 
and call it modified encounter model (m2). 

In order to distinguish between the two models, correlation times were measured 
at a constant temperature of 830 K and at different detector positions corresponding to 
different q and |g| for the single crystal and the polycrystalline sample respectively. 

MC simulations sites for different interaction energies between platinum atoms were 
carried out. One for no interaction between the atoms (si), one with strong repulsive 
potential between Pt atoms (s2) and one with an attractive potential between Pt atoms 
(s3) within the first neighbour shell (n = 1). 

The experimentally measured data as well as the data from MC simulations and the 
models for a diffusion mechanism in the single crystal are shown in figure [3] Note that 
the only adjustment parameter is the jump frequency, which is simply a multiplicative 
constant. The figure shows that the MC simulation with no interaction potential (si) 
complies with the encounter model for a system with no SRO (ml). The model with 
jump probabilities according to the encounter model and short-range order correction 
according to equation (|3]) (m2) gives the same values as the simulated data for strong 
Pt-Pt repulsion (s2) and describes the experimental data sufficiently well. For this model 
the fit produces a jump frequency of r^ 1 = 3.0(1) x 10 _3 s _1 . Unfortunately the quality 
of the measured data was insufficient to draw sound conclusions about solute-vacancy 
interaction. Obviously a model with only next-nearest-neighbour jumps (P 2 = 100 %, 
Pn^2 = 0) (m3) as well as a model with attractive force between Pt atoms (s3) is wrong. 

The polycrystalline sample was fitted by equation ([6]) and with the same jump 
probabilities as in the case of the single crystal. The short-range order correction 
was calculated by numerically integrating equation ((TJ over all orientations of q. As 
mentioned in section [3j this sample was only quasi-stable for a limited time. It can 
therefore not be assumed that a sufficient number of small crystallites was covered by 
the beam during the measurement. For averaging over a small number of crystallites 
with unknown orientation the Chudley-Elliott model [26J is, however, still the best guess. 
The calculated jump frequency using the modified encounter model as shown in figure 0] 
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(a) Radial scan of single crystal sample 
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(b) Azimuthal scan of single crystal sample 

Figure 3. Reciprocal correlation times in single crystal sample orientated in (110) 
direction at 830 K. The (001) direction lies at (j> = 0°. The radial scan (a) was taken 
at (j> = 33.2°. The azimuthal scan given in (b) was taken at 29 = 23° . Both datasets 
were fitted simultaneously with the classical encounter model with no short-range 
order (ml) and with the modified encounter model with correction for the short-range 
order calculated for strong Pt-Pt repulsion (m2). For comparison, a model with only 
next-nearest-neighbour jumps (m3) corrected for SRO is also shown. Monte Carlo 
simulations without interaction potential (si), with Pt-Pt nearest-neighbour repulsion 
(s2) and, again for comparison, with Pt-Pt attraction (s3) were carried out. 



was TQ~p C = 1.92(6) x lO -3 ^ -1 . The difference to the single crystal value arises most 
likely from a difference in sample temperature as a different furnace with a thermocouple 
at a different position with respect to the sample was used for this measurement. Again 
a model with only next-nearest neighbour jumps (m3) can obviously not appropriately 
describe the data measured. 

5.3. Activation energy and diffusion constant 

The activation energy is the sum of a vacancy formation and migration energy. The 
relationship between diffusivity, temperature and activation energy can be described 
fairly accurately by the Arrhenius law pQ. A number of measurements were performed 
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Figure 4. Reciprocal correlation times in polycrystalline sample at 830 K fitted with 
the classical encounter model (ml), modified encounter model (m2) and for comparison 
a model with only next- nearest- neighbour jumps (m3). 
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Figure 5. Chemical diffusion coefficients for polycrystalline (□) and single crystal 
sample (v) compared to results from tracer measurements [38] (O ). The green line 
(dashed) represents the fit for the single crystal sample, the blue line (solid) represents 
polycrystalline data and the red line (dash-dotted) represents the results of the tracer 
experiment. 





at different sample temperatures. Again measurements for different temperatures where 
repeated in order to guarantee thermal equilibrium of the sample. 

Using the Einstein relation |39j and the sum of partial diffusion coefficients into 
farther shells one can write: 




where |Aa n | is the radius of neighbour shell n. For fixed values of q one can use 
equation ([3]) to calculate the residence time To from the measured correlation times r(q) 
for each temperature. The lattice constant for Ni 97 Pt 3 was calculated using Vegard's 
law [40J and corrected for the thermal expansion at 830 K (a = 3.545 A). With this 
value and using the modified encounter model the chemical diffusion coefficients D can 



Studies of atomic diffusion in Ni-Pt solid solution by XPCS 



13 



be calculated for each temperature as shown in figure 

One should notice that a diffusion coefficient measured by the aXPCS technique is, 
owing to the coherency of the method, the chemical diffusion coefficient D. As already 
emphasized, exchanges between atoms of the same species do not influence the speckle 
pattern and are not recorded. However, in the case of our particular Ni-Pt solid-solution 
sample, chemical diffusion and self-diffusion coefficients of platinum are equal in a first 
approximation. This is simply a result of a low concentration of platinum atoms which 
makes direct exchanges between them highly improbable. 

Using the Arrhenius law the activation energy as well as the pre-exponential factor 
D can be calculated. Calculated activation energies for the single crystal and the 
polycrystalline sample, respectively were very similar with sc = 2.93(10) eV and 
Ea pc = 2.97(18) eV. These values are in excellent agreement with literature values for 
interdiffusion [3E], where the activation energy measured at much higher temperatures 
(between 900 and 1300°C) is (291.2 ± 3.7) kJ/mol = (3.02 ± 0.04) eV. The good 
agreement between our measurements and the data given in [38] allows us to conclude 
that there was no grain boundary diffusion contribution in the data given in that 
reference (apart from the point measured at the lowest temperature). The calculated 
pre-exponential factors are D sc = 1.5 (+4.6/-1.1) x 10~ 5 m 2 s _1 for the single crystal 
and D pc = 3.7 (+9.2/-3.4) x 10 _5 m 2 s _1 for the polycrystalline data. Again the value 
of Do = 9.2 (+3.4/-2.5) x 10 _5 m 2 s _1 found in [38] is in agreement with our results. 

6. Conclusions 

We could show that aXPCS allows the observation of low concentrations of solute atoms 
if the electronic form-factor contrast between solvent and solute atoms is high enough. 
On the other hand it is the question of the coherent synchrotron brilliance to be able 
to measure still more dilute samples. Correlated diffusion according to the encounter 
model with correction for short-range order effects is a simple atomic diffusion process 
and the measured factors are reliable and in very good agreement with the results 
of tracer methods. Diffusion rates measured by aXPCS are, however, in the order of 
10 -23 m 2 s _1 and therefore on the lower border of what can be reached by the best tracer 
measurements. For measuring diffusivities this low, aXPCS is the only method capable 
of determining atomic jump vectors. 

The Pt-Pt repulsion indicated by the Monte Carlo simulation is in agreement 
with the NiaPt phase [33J for low platinum concentration. Although for an atomic 
concentration as small as 3% one can not speak of a Ni 3 Pt phase, the characteristic 
atomic interactions which lead to its formation are present. It was, however, due 
to insufficient accuracy not yet possible to draw conclusions about solute-vacancy 
interaction from the comparison with simulations. With brighter synchrotron sources 
and better statistical significance coming along with them, it should be possible to 
investigate these features in the near future. Using Monte Carlo simulations on the 
bases of atomic interaction potentials and/or cluster expansion, aXPCS shows promise 
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in being a valid tool to confirm and test ab initio calculated potentials. 

X-ray photon correlation spectroscopy allows one to draw conclusions about the 
dominating diffusion mechanism. To distinguish subtle differences e.g. whether the 
encounter model or a model with slightly different jump probabilities takes place requires 
much better statistics than we achieved in our experiment. 

The experimental data from the polycrystalline sample showed generally good 
agreement in activation energy and diffusion rates with the tracer method. It seems also 
possible to reveal the dominating diffusion mechanism (in this case nearest-neighbour 
jumps) with these data. In order to prepare polycrystalline samples allowing for exact 
measurements of diffusion rates, efforts have to be made to prepare samples with 
sufficiently small crystallites and to eliminate crystallite growth effects. Especially 
aXPCS experiments on powder samples seem to be particularly promising in this regard. 

Overall we conclude that aXPCS is a very reliable method for investigating atomic 
diffusion for all types of alloys and provides a reasonable reference for ab initio 
simulations. 
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